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摘 要 : 准 稠密 匹配 是 多 视图 三 维 重建 的 重要 技术 ， 其 性 能 对 重建 结果 至 关 重 要 。 针 对 常用 的 Sif 算法 提取 的 种 子 点 
进行 准 稠密 匹配 正确 率 较 低 、 重 建 效果 不 佳 的 问题 ， 提 出 了 一 种 基于 尺度 不 变 Harris 角 点 特征 的 准 稠密 匹配 算法 。 该 
算法 首先 在 图 像 多 尺度 空间 构造 尺度 不 变 Harris 特征 ， 并 采用 余弦 距离 测度 对 不 同 视图 进行 双向 匹配 ; 然后 根据 稀疏 
匹配 获取 种 子 点 ， 采 用 最 优 最 先 匹 配 扩散 策略 进行 准 稠密 扩散 ; 最 后 采用 局 部 非 极 大 值 抑 制 策 略 对 匹配 结果 进行 
样 。 实 验 表明 ， 本 文 算法 提取 的 种 子 点 既 能 够 体现 场景 结构 信息 ， 又 具有 尺度 不 变 特性 ， 用 于 准 稠密 匹配 能 够 提高 
配 的 效果 和 精度 ， 是 一 种 有 效 的 用 于 三 维 重建 的 准 稠密 匹配 算法 。 
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Quasi-dense matching based on scale invariant Harris feature 
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Abstract: Quasi-dense matching is widely used in multi-view 3D reconstruction, and it is important for reconstruction results. 


和 Aiming at the quasi-dense matches diffused by the seed points extracted from Sift algorithm are less accurate, this paper proposed 


a quasi-dense matching algorithm based on scale invariant Harris corners. Firstly, it structured the scale invariant Harris features 


in multi-scale space, and the feature sets between different views are bidirectional matched by cosine distance similarity measure; 
Then the seeds selected from the initial matches are applied in quasi-dense matching algorithms by best and first propagation 
strategy; Finally, a local non-maximum suppression strategy is applied to resampling the quasi-dense matches. Experiments 
show that the seeds extracted by this algorithm can not only reflect the scene structure information, but also have scale invariant 
characteristics. And for quasi-dense diffusion, the matching effect and accuracy can be improved, and it is an effective quasi- 
dense matching algorithm for 3D reconstruction. 


Key Words: scale-invariant Harris feature; quasi-dense matching; local non-maximal suppression; 3D reconstruction. 


0 ”引言 准 稠密 匹配 算法 最 早 是 Lhuiller 等 人 "提出 ， 它 以 两 幅 图 

像 的 立体 稀 疏 匹配 作为 准 稠密 生长 的 种 子 点 ， 在 稀 玻 匹配 的 种 
随 着 计算 机 技术 的 飞速 发 展 ， 三 维 重建 技术 被 广泛 的 应 用 子 点 周围 进行 匹配 扩散 ， 得 到 稠密 匹配 。 因 此 稀 玻 种 子 点 获取 
到 影视 游戏 、 数 字 城 市 、 虚 拟 现实 、 机 器 人 导航 等 领域 。 基 于 对 准 稠密 匹配 和 重建 效果 至 关 重 要 。 近 年 来 ， 最 常用 的 稀疏 匹 
多 视图 的 三 维 重 建 因 其 重建 方法 便捷 及 适用 性 强 ， 成 为 三 维 重 配方 法 是 Lowe 等 人 BI 总结 的 SIFT(Scale-invariant feature 
建 领域 的 研究 重点 。 多 视图 三 维 重建 是 一 种 利用 同一 场景 的 多 transform) 算 法 , 由 于 SIFT 算法 具有 尺度 、 旋 转 和 光照 不 变性 ， 
张 不 同 视角 图 片 恢 复 三 维 模型 的 方法 ， 其 中 ， 立 体 匹 配 是 多 视 因此 在 初始 特征 提取 和 匹配 中 得 到 广泛 应 用 。 比 如 ， 张 需 等 人 
图 三 维 重建 的 关键 技术 。 多 视图 立体 匹配 主要 分 为 稀疏 匹配 、 中 利用 压缩 感知 理论 提出 了 一 种 基于 压缩 感知 的 Sift 准 稠密 匹 
眉 密 匹配 和 准 稠密 匹配 [1]， 由 于 准 稠密 匹配 计算 效率 高 而 且 能 配 算 法 。 王 伟 等 人 外 提出 了 一 种 基于 稀 琉 特征 点 的 匹配 扩散 方 
较 准 确 的 描述 三 维 场景 ， 近 年 来 一 直 是 研究 的 热点 。 法 用 来 估计 深度 。 赵 璐 瑞 等 人 四 通过 引入 彩色 信息 ， 提 出 了 一 
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种 基于 彩色 Sift 的 准 稠密 匹配 和 三 维 重建 方法 。 李 晓 明 等 人 区 
根据 结构 化 场景 特点 ， 提 出 一 种 基于 平面 单 应 约束 的 Sif 特征 
准 稠密 匹配 方法 。Guo 等 人 中 通过 改进 匹配 代价 函数 ， 提 出 了 
种 基于 Sif 的 准 稠密 匹配 方法 用 于 植物 的 三 维 点 云 获取 。 然 
而 ，S 这 算法 提取 的 特征 点 并 不 是 视觉 上 的 角 点 ， 特 征 点 在 图 
像 上 并 没有 实际 的 物理 特征 名， 并 且 耗 时 较 长 ， 会 影响 稀 玻 匹 
配 的 效果 。Harris 算法 四 提取 的 角 点 具有 较 好 的 信息 量 ,能够 很 


Ghinaxiy 合 作 期 乔 
孙 会 超 ， 等 : 基于 尺度 不 变 Harris 特征 的 准 稠密 匹配 算法 


然后 在 尺度 空间 中 检测 Harris 角 点 。 尺 度 空间 中 的 二 阶 和 矩 
阵 09 定 义 如 式 (C2) 所 示 。 


Livy,0p) Tily(x,y,op) 
Mnmonop)=abeg(lorj * a 


2 


RDGoyap) 有 Coyap) | 
其 中 : a 是 积分 尺度 ，op 是 微分 尺度 ，or=sop， 这 里 s 控制 a 
和 on 的 比值 , 在 这 里 取 s=0.8。 坟 和 工分 别 为 图 像 点 (x,y) 在 x 方 
向 和 y 方 向 的 方向 导数 ， 计 算 方向 导数 的 局 部 邻 域 是 原 图 像 与 


好 的 反映 物体 的 结构 和 轮廓 特征 目 ， 用 于 立体 匹配 会 获得 更 好 
的 效果 。 此 外 ， 对 于 多 视图 图 像 间 存在 的 尺度 问题 ， 
等 人 (9 基于 尺度 自动 选择 理论 ， 提出 了 Harris 尺度 不 变性 检测 
程 邦 胜 等 人 0 得 到 了 Harris 尺度 不 变性 检测 子 有 效 工作 的 
。 许 佳 佳 等 人 0 提出 了 一 种 结合 Harris 和 Sift 算 子 的 
图 像 配 准 算法 ， 提 高 了 图 像 配 准 精度 和 配 准 效率 。 
为 了 提高 准 稠密 匹配 和 三 维 重建 的 效果 ， 本 文 基于 
间 Harris 角 点 检测 和 Sift 特征 描述 方法 ， 提 出 了 一 种 适用 于 多 
视图 三 维 重 建 的 准 稠密 匹配 算法 ， 并 就 本 文 算法 和 最 有 具 代 表 性 
的 Sift 算法 进行 准 稠密 匹配 对 比 。 同 时 提出 非 极 大 值 抑制 重 采 
样 策略 对 准 稠密 匹配 结果 进行 优化 ， 以 去 除 元 余 。 本 文 整体 算 
法 框架 如 图 1 所 示 。 
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准 稠密 匹配 | | 


siHF 特 征 提取 与 稀疏 匹配 a 
输入 尺度 | [| [等 古 -| [最 机 准 得 密 
加 人 不 变 让 向 起 上》 向量 ES 区 | [和 > 三 维 
角 点 双向 生长 上 i 建 
检测 描述 匹配 匹配 扩散 重 
图 1 本 文 算法 框架 


1 SIHF 特征 检测 与 匹配 


Harris 算法 四 是 一 种 经 典 的 角 点 检测 算法 ， 在 没有 尺度 变 
化 的 情况 下 ，Harris 算法 对 图 像 的 旋转 变化 、 视 角 变 化 、 图 像 
噪声 具有 很 好 的 适应 性 [9。 为 克服 Harris 角 点 的 尺度 敏感 问题 ， 
本 节 引 入 尺度 空间 思想 外 检测 Harris 角 点 ， 并 采用 Sift 特征 描 


高 斯 核 函 数 o 卷 积 得 到 ,得 到 方向 导数 后 用 a 的 高 斯 核 函数 进 
行 平滑 。 计 算 Harris 角 点 响应 值 ， 本 文采 用 角 点 响应 函数 0 如 
式 (3) 所 示 。 


det(M (x, y,071,0D)) 
trace(M (x, y, 071,0p)) +eps 
其 中 :eps 是 一 个 极 小 量 (可 取 eps=10-)， 防 止 分 母 为 零 。 对 尺 
度 子 空间 的 每 一 幅 图 像 根据 式 (2) 和 G3) 计算 每 一 点 的 Harris 角 


G3) 


点 响应 值 ， 得 到 Harris 角 点 响应 金字 塔 ， 如 图 2 所 示 。 
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Harris 角 点 响应 金字 塔 


图 2 Harris 角 点 响应 金字 塔 


最 后 选取 尺度 不 变 兴趣 点 。 首 先 选 取 满足 阔 值 条 件 的 候选 
角 点 ， 当 Harris 角 点 响应 阐 值 设置 为 最 大 响应 值 Raax 的 
0.005~0.015 倍 , 角 点 检测 效果 最 她 ， 本 文 角 点 响应 阔 值 设置 
为 0.01Rmax。 然 后 ， 本 文 提出 尺度 空间 相 邻 层 圆 形 窗 非 极 大 值 
抑制 策略 提取 尺度 不 变 兴趣 点 。 如 图 3 所 示 ， 在 尺度 空间 中 ， 
在 候选 角 点 周围 选取 12 邻 域 的 圆 形 窗 ， 如 果 该 点 位 于 第 一 层 


述 方法 ， 构 造 出 尺度 不 变 Harris 特征 (Scale-Invariant Harris 
Feature，SIHF)， 并 将 SIHF 特征 用 于 稀 琉 匹配 。 本 节 有 具体 过 程 
如 下 。 

1.1 尺度 不 变性 Harris 角 点 检测 

首先 建立 图 像 的 高 斯 尺度 空间 。 将 原始 图 像 进行 下 采样 ， 
建立 一 个 图 像 金字 塔 ， 对 图 像 金 字 塔 的 每 一 层 通过 高 斯 滤波 建 
立 一 个 线性 尺度 子 空间 。 每 次 降 采 样 得 到 的 图 像 为 一 组 ， 每 组 
分 为 s 层 尺度 子 空间 ,本文 将 图 像 金字 塔 分 成 n 组 ,每 组 3 层 。 


对 于 一 张 分 辩 率 为 MxN 的 图 像 ，n 由 式 (1) 确 定 。 
n=logz{fmin(MxN)}-2 (]) 
对 尺度 子 空间 中 的 每 一 幅 图 像 ， 使 用 不 同 参 数 的 高 斯 核 进 


行 模糊 。 每 个 尺度 子 空间 包含 s=3 层 图 像 , 尺度 抽样 频率 为 ks= 
沪 ， 初 始 高 斯 模糊 因子 为 m=1.6。 原 图 像 为 高 斯 金字 塔 的 第 一 
组 , 第 一 组 图 像 的 高 斯 核 函数 为 a(1,s)=ooks*!, 高 斯 金字 塔 的 下 
核 函 数 与 上 一 组 的 中 间 层 相同 ， 尺 度 不 变性 
Harris 检测 子 的 微分 尺度 为 op=oln,s)。 


组 第 一 层 高 共 


和 最 后 一 层 , 则 在 相 邻 层 25 领域 进行 非 极 大 值 抑制 , 如 果 该 点 
位 于 的 中 间 层 , 则 在 相 邻 层 38 邻 域 进行 非 极 大 值 抑制 ， 当 满足 
上 述 条 件 时 ， 便 得 到 尺度 不 变 兴趣 点 。 


图 3 圆 形 窗 非 极 大 值 抑制 策略 


1.2 SIHF 特征 描述 和 双向 匹配 

得 到 兴趣 点 后 ， 对 兴趣 点 进行 特征 描述 ， 采 用 文献 后 的 方 
法 。 首 先 在 特征 点 周围 选取 16x16 的 邻 域 ， 计 算出 邻 域内 像素 
点 的 梯度 幅 值 和 梯度 方向 ， 并 统计 出 梯度 直方 图 ， 取 直方 图 峰 
值 方向 为 主 方向 ， 大 于 峰值 80% 的 为 辅 方向 ， 梯 度 幅 值 和 梯度 
方向 的 定义 如 式 (4) 所 示 : 
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mb)= Et Lb 2) + +d-Lony -1)? 
Ly+D) -Ly-) (4) 
L(x+1,y)—-L(x-1, y) 

其 中 : m(Xx, y) 和 0(x, y) 分 别 表示 梯度 幅 值 和 梯度 方向 。 然后 ， 
将 特征 邻 域 划分 成 4x4 的 小 栅 格 , 计算 每 个 小 栅 格 的 8 方向 梯 
度 幅 值 和 梯度 方向 , 这样 就 构成 了 128 维 的 特征 描述 子 。 最 后 ， 
为 了 去 除 光 照 变化 对 特征 点 的 影响 ， 对 特征 向 量 进行 归 一 化 处 
里 以 满足 光照 不 变性 。 
特征 匹配 常用 的 方法 是 欧 氏 距离 近邻 之 比 的 方法 。 由 于 特 
征 向 量 进行 了 归 一 化 ， 本 文采 用 工程 上 常用 的 余弦 距离 29 相似 
性 测度 进行 匹配 ， 余 弦 距 离 相 似 性 测度 计算 公式 如 下 : 


O(x, y) = tan 


下 
da; et 
sl (5) 


Di ; =arccos(d; ;) 


其 中 : xi、 久 分 别 为 两 幅 图 像 中 128 维 的 特征 点 集 ，d 为 两 个 
特征 向 量 的 余弦 角 , Dij 为 两 个 特征 向 量 的 余弦 距离 。 如 果 第 一 
个 特征 集中 的 特征 向 量 与 第 二 个 特征 集中 的 特征 向 量 最 近邻 
Dnearesi 与 次 近邻 D2nainearesi 的 比值 ratio 小 于 某 个 阔 值 ， 则 认为 
对 应 点 为 正确 匹配 。 为 得 到 更 高 的 匹配 正确 率 ， 采 用 双向 匹配 
03 策 略 ， 本 文 的 双向 匹配 策略 如 下 : 首先 ， 对 图 像 A 中 的 每 个 
特征 向 量 采用 穷尽 搜索 找到 图 像 B 中 满足 最 近邻 与 次 近邻 之 比 
的 特征 向 量 ， 匹 配 得 到 第 一 个 匹配 点 集 ; 然后 反 过 来 由 图 像 B 
到 图 像 A 匹配 得 到 第 二 个 点 集 最 后 ， 求 出 两 个 匹配 点 集 的 交 
集 即 为 最 终 匹 配 结果 。 


2 ”基于 SIHF 的 准 稠密 扩散 


准 稠密 匹配 使 用 稀 玻 匹配 作为 种 子 点 ， 再 在 种 子 点 周围 扩 
散 出 新 的 匹配 ， 得 到 稠密 匹配 。 本 文采 用 最 优 最 先 扩 散 策 略 争 
进行 匹配 扩散 ， 即 每 次 选取 一 个 最 优 的 种 子 点 ， 在 其 周围 进行 
匹配 扩散 。 然 后 根据 本 文 提 出 的 基于 局 部 非 极 大 值 抑制 策略 对 
准 稠密 匹配 进行 重 采样 ， 以 去 除 见 余 点 对 。 上 基体 过 程 如 下 。 
2.1 准 稠密 扩散 

得 到 稀 玻 匹配 后 , 首先 使 用 RANSAC 算法 加 获取 候选 种 子 
点 。RANSAC 算法 通过 重复 从 特征 点 数据 集中 获得 基本 子 集 ， 
利用 基本 子 集 估计 基本 矩阵 ， 得 到 满足 最 大 匹配 点 集 的 基本 拢 
和 车， 剔除 外 点 得 到 候选 种 子 点 ，RANSAC 算法 具体 过 程 参 考 文 
献 四 。 然 后 以 ZNCC(zero-mean normalized cross correlation) 得 分 
选取 种 子 点 ， 采 用 最 优 最 先 的 扩散 策略 ， 即 每 次 只 选取 一 个 得 
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视差 梯度 限制 的 候选 匹配 ， 并 计算 这 些 候 选 匹 配 的 ZNCC; 

0) 判 断 这 些 候选 匹配 的 ZNCC 分 数 和 置信 度 是 否 满足 阔 值 
条 件 ， 如果 这 些 候选 匹配 没有 被 匹配 过 (标志 位 为 0), 则 将 其 判 
为 正确 匹配 保存 下 来 (标志 位 置 1)， 并 放 入 种 子 队列 ; 

qd) 重复 以 上 三 个 步骤 ， 直 至 种 子 队列 为 空 。 

在 第 b) 步 中 ， 种 子 点 为、 尺 ' 的 (2N+1)x(2N+2) 邻 域 定义 如 
式 (6) 所 示 。 


N(x) = {uu—xe[l-N,N]} 


6 
N(x)={u,u -xe[-N,N]} (0) 
视差 梯度 限制 由 式 (7) 给 出 。 
N(x, Xx)={(u,u),ue N(x), ue N(x) (un) —(x-x) | < e} (7) 
其 中 : s 是 视差 梯度 限制 闪 值 ，(w, wi ) 是 扩散 得 到 的 候选 匹配 ， 


视差 梯度 限制 保证 种 子 点 对 到 候选 匹配 的 扩散 沿 同 一 方向 ， 选 
取 参 数 a=1， 邻 域 半 径 N=2。 


在 第 c) 步 中 ， 对 于 候选 匹配 (u,u') 采用 置信 度 限制 ， 置 信 
度 定义 如 式 (8) 所 示 。 
sC0)=maxflCc+A)-7TCOhAstf(GLO,(O+D]] (8) 


s"(w”) 的 计算 同上 式 , 置信 和 度 限制 就 表示 为 s(w) <t 或 者 s(x*) 
<t， 其 中 t 为 置信 度 阀 值 ， 防 止 匹配 扩散 到 灰 度 值 单一 ， 纹 理 
较 少 的 区 域 ， 一 般 取 {=0.01。 

2.2 局 部 非 极 大 值 抑制 

对 准 稠密 匹配 进一步 优化 不 仅 能 减少 元 余 点 集 数量 ， 优 化 
匹配 结果 ， 而 且 能 去 除 误 一 些 匹配 ， 提 高 点 集 集 质 量 。 本 文 提 
出 了 基于 ZNCC 得 分 的 局 部 非 极 大 值 抑制 重 采 样 方法 。 步 又 如 
下 : 


a) 建 立 与 原 图 像 大 小 相同 的 准 稠密 匹配 标志 矩阵， 用 来 存 
储 匹配 点 集 的 ZNCC 互相 关 分 数 ， 没 有 匹配 的 标志 位 置 0; 
b) 选 取 3x3 的 平滑 窗口 ， 判 断 每 一 个 扩散 得 到 的 匹配 点 是 
和 否 是 窗口 邻 域内 最 大 互相 关 点 , 若是 , 则 保留 下 来 , 否则 剔除 ; 
oO) 重复 以 上 过 程 ， 得 到 重 采样 后 的 匹配 点 集 。 通 过 对 准 稠 
密 匹 配点 集 重 采 样 ， 大 大 降低 了 点 集 的 数量 ， 去 除了 大 部 分 对 
重建 无 关 紧 要 的 元 余 点 ， 在 不 影响 重建 结果 的 情况 下 提高 了 重 
建 效 率 和 质量 。 


3 ”实验 结果 与 分 析 


本 文 实 验 环 境 为 Windows64 位 系统 , 运行 配置 为 Intel Core 
i5，2.5GHz CPU， 内 存 为 4GB， 算 法 采用 MATLAB R2017a 实 
现 。 实 验 包 括 两 个 部 分 ,分 别 为 SIHF 特征 稀疏 匹配 和 基于 SIHF 


分 最 高 的 种 子 点 进行 扩散 。 假 设 Cc, xi) 是 两 视图 图 像 中 的 对 应 
点 ， 将 匹配 按 ZNCC 分 数 从 高 到 低 排列 ， 称 为 种 子 队 列 ， 建 立 
两 个 与 图 像 大 小 相同 的 标志 矩阵 表示 对 应 像素 是 否 找 到 匹配 
(标志 为 1 表示 已 经 匹配 , 0 表示 未 被 匹配 )， 准 稠密 扩散 具体 过 
程 如 下 : 

a) 计 算 种 子 点 对 的 ZNCC 得 分 ， 将 得 分 从 高 到 低 排列 成 为 
种 子 队 列 ， 并 去 除 不 满足 阔 值 条 件 的 种 子 点 ; 

b) 从 种 子 队列 取出 得 分 最 高 的 种 子 点 ， 在 其 邻 域 搜索 满足 


特征 的 准 稠 密 匹配 。 
3.1 SIHF 特征 稀 玻 匹配 

为 验证 本 文 尺度 不 变 Harris 特征 在 稀疏 匹配 中 的 有 效 性 ， 
本 实验 选取 了 三 组 图 像 09 进 行 稀 玻 匹配 ， 并 与 Sift 算法 进行 比 
较 , 图像 之 间 变 换 关 系 已 知 可 以 客观 的 对 算法 精度 进行 评估 。 
图 4 为 三 组 原始 图 像 ， 分 别 为 尺度 变换 、 模 糊 变换 和 亮度 变换 
的 场景 。 首 先 对 图 像 分 别提 取 Sift 特征 和 SIHF 特征 ， 然 后 进 
行 特征 匹配 , 特征 向 量 的 匹配 阔 值 设置 为 ratio=0.8。 图 5 为 Sift 
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算法 的 匹配 结果 ， 从 左 到 右 依次 为 
变换 的 场景 。 图 6 为 本 文 SIHF 算 
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法 的 匹配 结果 ， 


次 为 尺度 变换 、 模 糊 变 换 和 亮度 变 ] 


(a) 尺 度 变 换 


换 的 场景 。 


(b) 模 糊 变 换 
4 原始 


(b) 模 糊 变 换 


尺度 变换 、 模 糊 变 换 和 亮度 。 果 均 为 原 


台 匹 配 ， 没 有 去 除外 点 。 两 种 方法 的 比较 数据 如 表 1 


从 无 到 右 依 。 ”所 示 。 
两 种 方法 匹配 结 


图 像 


(b) 模 糊 变 换 


5 Sift 算法 稀 琉 匹 配 


(b) 模 糊 变换 (b) 模 糊 变换 (b) 模 类 变换 
图 6 SIHF 算法 稀疏 匹配 
表 1 特征 提取 与 匹配 结果 
图 像 (a) (b) (9) 
方法 Sift 算法 ”SIHF 算法 ”S 进 算法 ”SIHF 算法 Si 进 算法 ”SIHF 算法 
总 匹配 对 数 2367 289 686 987 673 981 
正确 匹配 对 1608 253 277 926 434 924 
匹配 正确 率 ”67.93% 87.54% 40.38% 93.82% 64.48% 94.19% 
匹配 时 间 /s 。 359.56 112.62 334.09 130.17 157.78 60.09 
表 1 对 比 可 以 看 出 ， 对 于 一 般 场 景 的 图 像 匹 配 ， 本 文 算法 ”点 准 稠密 扩散 ， 得 到 87886 对 点 ; 图 8(b) 和 图 8 (c) 为 本 文 算法 
与 Sift 算法 同样 具有 很 好 的 尺度 适应 性 ， 而 且 计 算 效 率 也 得 到 ”” 准 稠密 扩散 ， 图 8(b) 为 原始 匹配 ， 得 到 89489 对 点 ， 图 8(c) 对 
提高 ， 本文 算法 与 Sift 算法 相 比 ， 对 不 同 模糊 和 复杂 光照 条 件 匹配 结果 进行 了 非 极 大 值 抑制 重 采 样 (以 下 简称 重 采 样 )， 得 到 
的 彩色 图 像 能 够 提取 出 更 多 稳定 的 角 点 特征 ;在 特征 匹配 上 ， 17208 对 点 。 
本 文采 用 的 余弦 距离 度量 的 双向 匹配 策略 ， 进 一 步 剔 除 大 部 分 为 了 更 直观 地 体现 本 文 准 稠密 匹配 算法 的 优越 性 ， 并 从 客 
的 误 匹 配点 和 重复 匹配 点 对 ， 能 够 得 到 更 高 的 匹配 准确 率 。 观 上 评价 算法 性 能 ， 将 准 稠密 匹配 进行 两 视图 三 维 重建 ， 采 用 
3.2 SIHF 特征 准 稠密 扩散 重建 后 的 三 维 点 云 在 图 像 上 的 重 投影 误差 做 度量 准则 。 如 果 重 
为 验证 本 文 准 稠密 匹配 算法 的 有 效 性 ， 实 验 采 用 普通 相机 拍摄 ”投影 误差 满足 阔 值 条 件 ， 判 断 匹 配 是 正确 的 ， 重 投影 误差 定义 
的 一 组 Box 场景 ， 图 7(a) 为 600x800 的 原始 图 像 。 实 验 首先 对 ”如 下 : 
图 像 进行 稀疏 匹配 ， 然 后 进行 准 稠密 扩散 ， 稀 玻 匹 配 中 特征 向 有 
量 的 匹配 阔 值 设置 为 ratio-0.8。 图 7(b) 为 Sift 算法 稀疏 匹配 结 CN ® 
果 , 得 到 265 对 点 ; 图 7(c) 为 本 文 SIHF 特征 稀 疏 匹配 结果 , 得 。 式 (9) 中 ，wW 和 ws 分别 为 重建 后 的 三 维 点 根据 投影 关系 映射 回 
到 295 对 点 。 图 8 为 准 稠密 扩散 结果 ， 其 中 图 8(a) 为 S 进 种 子 图像 平面 的 坐标 ， 码 和 轻 为 匹配 点 实际 的 图 像 坐标 ，drepgror 表 
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录用 稿 和 孙 会 超 ， 等 : 基于 尺度 不 变 Harris 特征 的 准 稠密 匹 

示 重 投影 误差 。 数量 大 大 提高 。 其 中 图 9%(@) 的 准 稠密 重建 结果 经 过 了 局 部 非 极 
准 稠密 重建 采用 经 典 的 SFM 算法 073， 准 稠密 重建 的 结果 ”大 值 抑制 重 采样 ， 得 到 点 云 个 数 为 2886， 可 以 看 出 经 过 重 采 样 

如 图 (8) 所 示 。 图 9(a) 是 常用 的 基于 Sif 特征 的 准 稠密 重建 结果 ， 的 点 云 保留 了 模型 关键 位 置 的 点 ， 去 除了 大 部 分 对 建 模 无 关 紧 

得 到 三 维 点 云 个 数 为 2993; 图 9(b) 是 本 文 算法 准 稠密 重建 结果 ， 要 的 点 。 两 种 方法 的 具体 比较 数据 如 表 2 所 示 。 

得 到 三 维 点 云 个 数 为 13695， 可 以 看 出 本 文 算法 的 有 效 重建 点 


(a) 原 始 图 像 


(a) Sift 


(a) Sift 


表 2 算法 性 能 比较 


图 9 准 稠密 


(c) SIHF+ 重 采样 


(b) Sift 
图 7 稀 玖 匹配 


(b) SIHF 
图 8 准 稠密 匹配 


(c) SIHF+ 重 采样 


8; 强 吕 


(b) SIHF 
三 维 重 建 


(c) SIHF+ 重 采样 


影响 ,提出 了 一 种 基于 尺度 不 变 Harris 特征 的 i 


侍 稠 密 匹 配 算法 。 


在 稀 玻 匹配 过 程 中 ,本 文 基 于 Harris 尺度 不 变性 角 点 检测 算法 ， 


稀疏 。 ”” 准 笛 密 重建 本 
2 4 构造 出 重 棒 的 特征 描述 子 用 于 稀 玻 匹配 ， 实 验 表明 本 文 的 稀疏 
i 0 匹配 算法 具有 较 好 的 效果 。 在 准 稠密 匹配 过 程 中 ， 提 出 了 基于 
有 WE SIHF 特征 的 准 稠密 扩散 方法 , 同时 提出 局 部 非 极 大 值 抑制 重 采 

RE 样 策略 优化 匹配 结果 。 实 验 表明 ， 本 文 的 准 稠密 匹配 算法 能 
大 幅 提 高 准 稠密 匹配 正确 率 ， 用 于 三 维 重建 能 够 获得 较为 满意 


从 表 2 可 以 看 出 , 基于 SIFT 特征 和 基于 SIHF 特征 


密 扩 散 都 能 扩散 到 全 局 有 纹理 


的 区 域 ， 但 是 重建 结果 经 过 重 投 。 得 到 较 好 的 准 稠 密 匹配 ， 下 一 # 


的 结果 。 本 文 算法 在 准 稠密 扩散 过 程 中 对 于 复杂 纹理 区 域 能 够 
将 重点 研究 如 何 提高 简单 纹理 


的 准 笛 


影 去 除 误差 较 大 的 匹配 点 后 ,本 文 算法 得 到 的 匹配 正确 率 更 高 ， 
重建 的 点 云 质 量 更 好 。 从 图 9 Box 场景 的 三 维 点 云 可 以 看 


三 维 


区 域 的 匹配 精度 和 匹配 数量 。 
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